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РЕШЕНИЄЕ 3р ЗАДАЧИ КОМПЬЮТЕРНОЙ ТОМОГРАФИЙ ПО 
ИЗВЕСТНЬЮМ ТОМОГРАММАМ НА СИСТЕМЕ 
ПРОИЗВОЛЬНЬІХ ПЛОСКОСТЕЙ 


У статті будується та досліджується метод відновлення внутрішньої структури тривимірного тіла 
за допомогою поліноміальної інтерфлетації з використанням відомих томограм (слідів), що лежать 
на системі довільних площин, який є узагальненням методу відновлення тіла за відомими 
томограмами на системі трьох груп паралельних площин. Сформульовані та доведені теореми про 
інтерфлетаційні властивості та похибку побудованого оператора. 

Ключові слова: комп'ютерна томографія, інтерфлетація, відновлення. 


Тре агісіє 15 Базед гесоуегу теїпод апа іпуезбісатед Ше «ргее-Фітепзіопа! іпіета! 5ігисіцте ої "ре Боду ц5іпе 
роїупогаіа! іпіегйатайоп цзіпо Кпомуп іоплодгатя (асе5) Іушпе оп а 5узіега ої гапаот ріапе5, мрісі і5 а 
еепегайганоп ої Фе плеїпод ої гезіогіпо, Ше Боду Кпомуп іотовгатя оп а 5узіеті ої геє отгопр5 ої рагаїсі 
ріапез. Копииіата апа ргоуса пеогетя оп іпіегЙеіаїсіупі ргорегіїез апа бий орегаїог еттог. 

Кеу ууогаз: согариїег готортарбу, іпіегПагайоп, гезіогайоп. 


В статье стройится и исследуєтся метод восстановления внутренней структурьшю трехмерного тела с 
помощью полиномиальной интерфлетации с использованиєм известньх томограмм (следов), лежащих 
на системе произвольньвх плоскостей, которьй являєтся обобщениєм метода восстановления тела по 
известньтм томограммам на системе трех групп параллельньх плоскості ей. Сформулированьт и доказаньт 
теоремьт о интерфлетационнькх свойствах и погрешностях построенного оператора. 

Ключевьге слова: компьютерная томография, интерфлетация, восстановлениє. 


Вступ 


В останній чверті 20-го століття в практиці медичних досліджень, а також 
при неруйнівному контролі тривимірних об'єктів, при проведенні наукових 
досліджень у різних галузях науки 1 техніки тощо, знайшли широке застосування 
комп'ютерні томографи (1-31, які дозволяють відновлювати внутрішню структуру 
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тіла не розрізаючи його. При цьому виник новий клас задач - задач відновлення 
внутрішньої структури тривимірного тіла за відомими його томограмами на 
декількох площинах. 

У роботах |4-5| були побудовані та досліджені оператори поліноміальної та 
сплайн-інтерфлетації 16) функції трьох змінних за відомими слідами на системі 
трьох груп площин (у кожній групі площини паралельні), 1 на основі цих 
операторів була розв'язана задача тривимірної комп'ютерної томографії, у 
випадку, коли відомі томограми в системі трьох груп перерізаних площин 
(паралельна схема сканування). Була доведена висока точність розроблених 
методів. Відомими є методи розв'язання 3) задачі комп'ютерної томографії за 
допомогою конусоподібної та спіральної схеми сканування, в яких задані 
площини не є паралельними |7-9. 

Робота присвячена розв'язанню задачі відновлення внутрішньої структури 
тривимірного тіла за відомими томограмами, що лежать на системі довільно 
розташованих площин. Розроблений в статті метод є узагальненням методів, 
розроблених авторами в роботах |4-5| та є більш точним, ніж відомі методи 
авторів |7-91. 


Побудова оператора інтерфлетації на системі довільно розміщених 
площин 


Нехай скалярнозначна функція | /(х), х з (х,хо,хз), хє В) описує деяку 
фізичну характеристику внутрішньої структури (наприклад, щільність, коефіцієнт 
поглинання тощо). Джерелом інформації про функцію  /(х), тобто про 
внутрішню структуру тривимірного тіла, будемо вважати набір будь - яких М 
перерізаних площин, а також набір томограм Т, (х) , які лежать на цих площинах, 


які задаються наступними рівняннями 


3 
По (х)я а -ураб0,Каі, 
ре 
Вважаємо, що в одній точці перетинається на більше трьох площини. 
Введемо наступні позначення: 


о 
1. 0 та З|ад 4 | - вектор, направлений вздовж лінії перетину 
дро Фр з 
площин 05, -0, 0, -0);; 
2. Мік, І) п, СП,ОПу ЗР З (ХульХочо Хз) 9 Фіто, | де Р 


точка перетину трьох площин; М - множина точок перетину; 
3. Ги-П ОП, яд, - ребра, за якими перетинаються дві площини, на яких 


лежать відповідні томограми; 
бар біо б 

4.00 ДЗ 4 4з|- вИЗНачник, складений із коефіцієнтів рівняння 
Чі 4 Є 


заданих площин; 
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5. Т. (х)-томограма, яка лежить на площині П,. 


Визначення 1 |4)|. Томограмою т.к) (слідом функції /(х) на площині 


о(х)-0 за умови, що коефіцієнти а,,і «1,3 не дорівнюють нулю, будемо 
називати одну з трьох функцій 
Лого) | ЛО оо з / ар Хо о)» у 0 
Т(х) 24 (ооо) ЗО (У паху а) / оз з), ар 0, 
ЛО ооо (ьо) Ло (У р ар Хо) / дз) аз 0 
(хх), ха -0 
хач(о,хз), хо є0 
(х»,Хз), хі 0 
де Ху Хр» -13 - вирази, що отримуються розв'язанням рівняння й), (х) - 0 
відносно змінної Хр: 


Нехай томограми Т,(х), Т, (Х), Т, (хХ) перетинаються в точці Й,,. Позначимо 


к Ті Ти Ти 
и 92 Рувото- Ф(х) ач (х), у (х) є Доу о (ХУ). 
жі ні ні 


Теорема 1. Для існування функції /., (х) є С"(О) із заданими томограмами 
т 2)із І, п, т, (2), Ка І, т, то (х), Ге 1,5, для якої виконуються умови 
Ім СОЇ, З ТАС о ОЇ, З РОЇ, 0 Фк ОЇ У ТиСО, з 
необхідно та достатньо, щоб сліди Т, |(х), 4 -1,2,3, Фі, із іп, Каеїт, 
І -1,5 задовольняли умові Т. 4(Х) є С"(В?), к20 та умовам С.М, Нікольського, 
які на ребрі Г,, зводяться до перевірки рівностей 
Т заці (х)) -Т о (ші (2) 


тобто, значення томограм на лінії перетину повинні співпадати для всіх томограм, 
що перетинаються. Аналогічний вигляд мають ці умови на ребрах Г, Г,,. 


о3(юе0 од дао? 


У точці /,; умови Нікольського зводяться до перевірки рівностей 
І І Е 
Т ушу (2) -Т (ух) - То) 


тобто, значення томограм у точці перетину повинні співпадати для всіх томограм, 
що перетинаються. 
Доведення. Оператор 1 (х) побудуємо у вигляді 


І (Х) З по ) ха ря кій у сх па ри Ре ої - їі мб | Ста ), х, 0) 
Іо Та З Ту и (2) є (цу Сх), 


ої; (х)-0,2, (хує0 3 (хе, (х)-0 од (буси (хує0 
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иа ро ОО о р нд в-23 4 і В 
Аналогічно визначаються оператори Риб Нофо, Па ; 
У роботі С.М. Нікольського |10) доведено, що для того, щоб існувала функція 
Лого) є С" (0), яка має сліди 
Р; (Ха) 


ся фу 823, рабу, 


хт0 
необхідно та достатньо, щоб граничні функції ф,. і (2) задовольняли в точці (0,0,0) 


умовам вигляду 


д?ф ах 9"ф, зх 0,г 
дф. а 2 09,80) ; 5; р- 12,3, аВеб0,г 
дхі ху 
Хрт0 й хуя 


Цей факт дозволяє будувати оператори /., (х) у вигляді (2). 
Теорема 1 доведена. 
Розповсюдження результатів на випадок більшої кількості площин ніж три 
проводиться в теоремі 4, завдяки використанню розкладу одиниці. 
Теорема 2. Нехай внутрішня структура тривимірного тіла описується 


функцією  /(х) є С"(0) (г 23), яка має томограми Т,(х), Кк -1,МУ, задані на 


площинах П, відповідно, та задовольняє умови / сода б 


020). Тоді для 
похибки К.,, / (х) - (1 - 2.) /(х) наближеного відновлення внутрішньої структури 


ЛО) оператором /.,(х), побудованим за допомогою даного набору площин та 


томограм, виконується Уч 


і; 2, 93; 
т 
Ки / (9 з ї І Парора г б 77 пС ор фр в 17 Фа. (3) 


Ді Пк їКІ 


Доведення. носу до уваги тотожність 


(1-70) 2(І- їі -Пи -ї т ГАЮ і; ЛНВ с дї т ВР ОГО з 


-(-вХі-вуі-щуоо 


а також тотожність типу 


5 () 2,93 3,1 і, 2 
дЇ Ти іє Тік 
ін Ж ус ви й. 

Кі 7 їКІ 


ща Фі. 
пі о З Чу : : К 
ідставимо цей інтегральний вираз, а також аналогічні вирази для (І -Ї ) Д(), 


(І - ди ) 702) в попередню формулу 


Ки (9 - (1-Їа )(9 та (Г- 0 Хід Х-дя ус за 
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й 23 т92,03 то3ої л9о2 
-РМЦІ| , Л Пи є З ен Ні, | 41, Ай, А). 
о 0 0 


ні ПК їКІ 


Теорема 2 доведена. 


Теорема 3. Нехай множина довільних томограм, які знаходяться на 
площинах, що задаються рівняннями: 


п, :о0 (2) 20,к -1,М, 


задовольняє умову: в одній точці Й,, - ПІ, о П2, г ПЗ, перетинаються не більше 


І. 
трьох томограм; 
Тоді система функцій 


М 
П оо 
де 

ян 0 з дак 


М 
П и) 
й 
дні ЮЇ 

має властивості | й (Йурр) Є б; удррду» із -І,п, К,К'зІ,т, 1, «15. та | є 


т 
розкладом одиниці, тобто Я ли (2)21. 


і, ЄМ 
іяКа 


Доведення.  Врахуємо, що У/«-у-ПІі,ОП2,СП3, тобто  зі(Й/ и. )-0, 
о2(Й) 2 0, а/д -0. 


П сої (Й) П 2) П 3 у) 
уч уні уч 


оо ре укі ек укі Чо із ке, аї Ми 
ПРООНКФУО І сот о ; ч 
щ Щ 5 0, в інших випадках 
Поу) П од (у) Пози) 
узі узі уч! 
зі жк ні 


абидуиді» пд , юка, 1,1 213. 


Теорема 3 доведена. 


Лема 1. Нехай  томограми Т,, (Х), і -1,п, т, Кк -1т, То (Х), / -1,5 
задовольняють умовам теореми 2.9, тоді система функцій Л, (х) є розкладом 
одиниці: 


5 


т п 

2, 3, Ху идз! 
ї-і ло ПЕ 

їн Ї о Кі 0 Ізі К 


Доведення. Рівність для й, (х) можна записати у вигляді 
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По, 9. Под», 9 Деиз 9 
ці ю Вч 


їм (х)з з ужі ук Ве - 
Поїх, (Ру) Под, (Й) Под (7, й 
Ра і б. 


по пі, (9) ут спо, (Х) ру пз, 9) 


нат Фіти, ( Й) у з оди (Й) 2 Зп (Й) 
НЕ 


Тепер просумуємо цей вираз 


п т 5 

2 УЗ З йн 9) є 
ї«Р (Ко І-І 

іяКаї Кізі ІзКі 


У 3 У й іт, (9 Й оп, (9 Й оЗпз | 


ТО КА ІЗ неї іі, (Йи) ул оп», (и) Ва Зп, ( Ра 
іні Каізе! Гжіжк| неі ук Ві 
п Фіг, 0) о3пз, (9 


-У|П Р ЛЬ Ре 
їзі | цеї іти, (Ри) | єторул оп», (Ри) 8-1 3 (Йи) 
іяк|| ужі Кізі | уж Гекскі | Ве! 


Тут враховується, що 
по Фі, (0) ру оц (х) у Зла, 9) 


узі Фіг, (Й ор) У ус з оп», (Й ПО і 1 1 Зп (7, ВУ 
ижі 


є базисними поліномами Лагранжа. З цього виходить, що 


п.|тооіц, (х) т | то суда» (Х З пз (Х) 

5 Пан 3 Пн Пн 
їз | неї іп, (Йи) ке | узі оп», (и) і-ї | Ве! Зп (Йи) 

іяісеЇ| нжі Кеті | ук Гекті| Ве! 


ПІ 


Лема доведена. 


Теорема 4. Нехай томограми Т, (х) є С" (К/), р23 задовольняють умовам 
С.М. Нікольського на ребрах і в точці перетину площин. Тоді функція 


ІЩх)а 23 й СД (0) 


(уЄТем 
є поліноміальним інтерфлетантом із властивостями 
Щх) є СО), ЩО, - ТО), 8-1, (4) 


При цьому У/(х) є С"(О), г 23, що задовольняє умовам теореми 3, виконується 
рівність 


Щху ЗО), Го ОВО), ВОГО У ОВУ ГОЮ 


(Гм 
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де К.//(х) визначається формулою (3), а й, (х) - допоміжні поліноми, що 


визначаються в теоремі 3. 
Доведення. Згідно з формулою (2) при побудові операторів /., (х) 


використовуються функції 7, (2), Т»., (2), Ту (х) є С" (7), р23. 


Тому Щх) є С"(К?) або, якщо СО. с 87, то Щ(х) є С"(О), Доведення співвідношень 
(4) проводиться так (використовуються умови теореми 1 і те, що система функцій 
й (9) є розкладом одиниці): 


огід «У У У сдіцисої | о У УмооЇ, ІОдЇ, ТОЇ 


їзІ Каї ІчІ п; ізі Каї І-І 


Враховуючи, що У У У Ам) - І, формула (дви ГО9 2 709, 


ізі Ка І-і 
ік о Кін  Інкоі 


перепишеться наступним чином 


Лог У У Упи (2) ДРУДЕ; б Ви СОЇ 


ізо ке А 
Звідси отримуємо 


п 5 


то У У УпиобінГОО КУ У У иодву Од з (000 Вод) ГО). 


їз І Ка І-І їз К-і І-і 
Теорема 4 доведена. 


Теорема 5. Абсолютна неусувна похибка Е побудованого інтерфлетанта в 
припущені, що /(х, у,2) на площинах П,, тобто відповідні томограми задані 


наближено 6,, тобто 


то-Т ОЇ є бу ел, 


а також 


Т оо, - РО9Ї, 


«бі» Каїп,ізіт, 


ін Сх) -РО) 


о (х)є0, у (хе 0 


«би» 


о 2020, 0 (2)20 
дорівнює 
Ес У букбунд кб, нд кб дб 
(бю Пєм 
Тестовий приклад відновлення функції трьох змінних за допомогою 
оператора інтерфлетації за відомими її слідами на системі довільно 
розташованих площин 


За викладеною методикою був розроблений комплекс програм у системі 
комп'ютерної математики МаїйСаа. Результати її тестування демонструють 
високу точність. 

Продемонструємо результати роботи програми. 

Нехай задані 4 площини П, : 05, (х) -0 у вигляді 

лявумачті сей. саке с ноу - 4 


ро оо У о о о 
1 В 2 58 3 В 4 «З 
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які перетинаються в чотирьох точках: 
Й зз 7 (0,0,1), Й а 2 (0,-1.5,2.5), Йза 2 (-1.5,0,2.5), Роза 2 (-1.5, -1.5,1). 
Як бачимо, в одній точці не перетинається більше трьох площин. 
Про функцію К()) є ху я 10х, -5ху відомі її сліди на заданих площинах, 


тобто томограми 


о) щ ВО зхі-б5хж5х 5, 
т.е ко, аю хонбх НІ15Х;кУ; 
т, (хз раге ЗНЯЛИ оно УЮ 
то)з КО о хі-нбх чн 15х, 45. 


За теоремою 4 був побудований оператор інтерфлетації, який використовує 
лише задані томограми та рівняння площин, на яких лежать томограми. Після 
спрощення побудований оператор набуває вигляду 


Годо хм 
Висновок 


З тестового прикладу робимо висновок, що побудований оператор 
інтерфлетації за відомими  томограмами (слідами) на системі довільно 
розташованих площин точно відновив квадратичну функцію, чого не можливо 
досягнути за допомогою операторів інтерполяції, які використовуються в 
сучасних методах комп'ютерної томографії. 

Викладений метод відновлення внутрішньої структури тривимірного тіла за 
допомогою операторів інтерфлетації у випадку відомих томограм, що лежать на 
системі будь-яких площин, є узагальненням розробленого авторами методу 
відновлення за відомими томограмами на системі трьох груп перерізаних площин 
і має таку ж високу точність. 
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КО не 


оо 


ВЕ5СМЕ 


У.І. Рег5пупа,О.У. 51ййп,У. О.Разісппік 
5о0Їуїпо Ше ргобіет ої 3р сотриїег готогсгарпу Їог пе Кпоууп апа ії5 5у5кет ої 
агрікгагу тотоєгатя оп ріапе5 

Тіз агіїсіе 15 сопзігисіед апа іпуе5іїєаїедй таїбетайса! піоде! ої (ууо- 
дітепз5іопа! согариіед їотовбгарпу и5іпє іпіегбаїайоп Іипсбоп5 об (геє уагіабіез 
Ккпомтп іотогтат5 5їидеа Боду. І 15 Бепеуеай Фаї іотоєгатя аге бїуеп оп а 5у5ієт ої 
агбіїгагу ріапез (аг опе роїпі іпіег5есіє аї пло5і їргее ріапе5). ТПе агіїсіе рге5епіє Ше 
сопсері ої (отоєтарбу їп Бе птааїпетайса! 5еп5е, а5 Ше їгасе об а Типсіїоп ої Фгее 
уагіабіе5 їп Бе бїуеп ріапе апа КФРе аїсогіїрта ої ігапзіайоп котоєгарпу ітаєе а5 а 
Дипсбоп, Фе агештепіз ої уУрбісп 15 Ше пигбег ої ре Поиге апа Ше соогаїпаіез ої 
ріхеї8. Тіз біуе5 Фе орропипіїу о м/огк міїп (оплоргатзя аз Гапсйопе, 1.е аПом'пе ог 
Нбосшге пштбег іо обіаїп її5 плаєбе апа рБієвПпеБі соїог сотропепі аї Фе з5ресійед роїпі 
ої дгаміпо. 

М/е сопзігисі апа 5їиду Ше орегаїог іпіегПагайоп Гипсіоп ої Шгее уагіабіе8 
Ккпомт Ве БоПомтпеє Гапсйїоп5 оп а 5у5іет ої агріїгагу ріапез. Ме ргоме а Шеогет оп 
Бе депега Бог ої еггог арргохітайоп ої а Гапсіїоп ої ігее уагіабіе5 соп5ігискей 
іпкегбакайоп орегаїог іп іпіеєгаї Їогп. АЇ5о біме ап е5іітаїе шпгесоуегабіє еттог 
ехрегітепіа! Чака. Депоп5ігаїсд а (е5і са5е Їог Ме соп5ігисіїоп ої Ше орегаїог 
іпгегйакайоп дпадгайс Гипсйоп, апа її Ба5 Бееп 5помитп Шаг їбе соп5ігисіва 5агисіите 
Ьгіпо5 ехасіу Віз Гипспоп, урісП 15 пої ре сіа8віс іпіегроїайноп орегаїог8. 

Тре ргорозед плеїбой еззепіаПу дїйег5 йот Ше ехі5ійиє опе5, їБаг її сап Бе 
сагтіед оці ргосе8зіпе готобтатз5 Шаї до пої Пе їп рагаїеі ріапез (е.є., п Фе зітріезі 
са5е, Ше їоторбтат тау Бе агтапеей 5узіет ої ргее єгоипр5 ої ріапез5 рагайе! о Ше 
соогдїпаїє ріапе8). ТПе плеїбдод піакез їс роз55ібіе 10 50Їуе Ше ргобіет ої іргее- 
Фітепзіопа! согариіедй їотоєгарпу їог а Жіпдатепіайу пеуу 5сПетеє Їог аа 
соЙеспоп. Ког ехатріе, її регтпії5 Ше ц5е ої Гап-соПесіоп 5спетез ої іпЇогтайоп їп 
еаср ої Фе ріапез їп убісі Пе їБе котоєбтат. 
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Ю.І.Першина,О.В. Шилін,В.А. Пасічник 
Розв'язання ЗР задачі комп'ютерної томографії за відомими томограмами на 
системі довільних площин 

У даній статті побудовано і досліджено математичну модель двовимірної 
комп'ютерної томографії з використанням інтерфлетації функцій трьох змінних за 
відомими томограмами досліджуваного тіла. Вважається, що томограми лежать 
на системі заданих довільних площин (в одній точці перетинається не більше 
трьох площин). У статті дається поняття томограми в математичному сенсі як слід 
від функції трьох змінних на заданій площині і побудований алгоритм переводу 
зображення томограми у функціональну залежність, аргументами якої є номер 
малюнка і координати пікселів. Це дає можливість працювати з томограмами, як з 
функціями, тобто дозволяє за номером малюнка отримувати його зображення 1 
виділяти компоненту кольору в зазначеній точці малюнка. 

У роботі будується і досліджується оператор інтерфлетації функції трьох 
змінних за відомими слідами функції на системі довільних площин. Доводиться 
теорема про загальний вигляд похибки наближення функції трьох змінних, 
побудованих оператором інтерфлетації в інтегральному вигляді. Також приведена 
оцінка неусувної похибки експериментальних даних. Продемонстровано тестовий 
приклад побудови оператора інтерфлетації для квадратичної функції, і показано, 
що побудована конструкція наближає цю функцію точно, чого не можна сказати 
про класичні оператори інтерполяції. 

Запропонований метод істотно відрізняється від існуючих тим, що в ньому 
може проводитися обробка томограм, що не лежать в паралельних площинах 
(наприклад, у найпростішому випадку томограми можуть розташовуватися у 
системі трьох груп площин, паралельних координатним площинам). Метод дає 
можливість вирішувати тривимірну задачу комп'ютерної томографії для 
принципово нової схеми збору даних. Наприклад, він допускає використання 
віяльної схеми збору інформації в кожній з площин, в яких лежать томограми. 

Надійшла до редакції 01.09.2015 
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